Electronic Transport Through EuO Spin Filter Tunnel Junctions 
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Epitaxial spin filter tunnel junctions based on the ferromagnetic semiconductor europium monox- 
ide, EuO, are investigated by means of density functional theory. In particular, we focus on the spin 
transport properties of Cu(100)/EuO(100)/Cu(100) junctions. The dependence of the transmission 
coefficient and the current-voltage curves on the interface spacing and on the EuO thickness is ex- 
plained in terms of the EuO density of states and the complex band structure. Furthermore we also 
discuss the relation between the spin transport properties and the Cu-EuO interface geometry. The 
level alignment of the junction is sensitively affected by the interface spacing, since this determines 
the charge transfer between EuO and the Cu electrodes. Our calculations indicate that EuO epi- 
taxially grown on Cu can act as a perfect spin filter, with a spin polarization of the current close to 
100%, and with both the Eu-5d conduction band and the Eu-4/ valence band states contributing 
to the coherent transport. For epitaxial EuO on Cu a symmetry filtering is observed, with the 
Ai states dominating the transmission. This leads to a transport gap larger than the fundamental 
EuO band gap. Importantly the high spin polarization of the current is preserved up to large bias 
voltages. 

PACS numbers: 71.20.-b, 71.20.Be, 75.50.Pp, 72.25 Dc 
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I. INTRODUCTION 

Magnetic tunnel junctions using MgO as a tunnel bar- 
rier generally display large tunneling magneto resistance 
(TMR) and they are now widely used in read heads for 
ultra high-density hard-disk drives as well as in random 
access memory devices^. Currently a room temperature 
TMR ratio as high as 500% can be reached in junctions 
made of polycrystalline Fe (more precisely CoFeB) sand- 
wiching a thin MgO insulating layer—. Despite such huge 
TMR ratios the Fe/MgO system presents two limiting 
aspects for future applications. Firstly, there is no flexi- 
bility on the materials side, which means that high TMR 
ratios are achieved only for the Fe/MgO combination and 
only for a particular crystal orientation, (100). Further- 
more high crystallinity is needed as the electrodes them- 
selves carry only a moderate spin-polarization of the con- 
duction electrons. Secondly, the TMR consistently de- 
creases with increasing the applied voltage and/or the 
temperature and so the operation of a practical device 
utilizing a magnetic tunnel junction is limited to low 
bias 3 -. 

The origin of the high TMR found in Fe/MgO-based 
junctions is rooted in the spin filtering effect^ that MgO 
exerts on the conduction electrons injected from the Fc 
electrodes^—. Only electrons with specific symmetry 
and with small transverse momentum contribute signif- 
icantly to the current. For energies around the Fermi 
energy, Ep, these are found only in the majority Fe 
states. As a consequence of this symmetry filtering the 
decay rate of the conduction electron wave-function into 
the MgO barrier of the majority spins is smaller than 
that of the minority ones^. These decay rates are de- 
termined by the complex band structure (CBS) of the 



tunnel barrier^. A crucial feature of the symmetry spin 
filtering effect is that TMR ratios larger than 1000% can 
be achieved, despite the fact that the spin-polarization 
of the electrodes is only around 65% (for CoFeB). 

Another possible strategy for obtaining a spin filter ef- 
fect, and thus very large TMR ratios, is that of using 
a ferromagnetic insulator as spacer between the metallic 
electrodes. In this case the energy barrier has a different 
height for majority and minority spins, which leads to 
the suppression of the current for one of the two spin- 
species for thick enough junctions. The tunnel current 
density is proportional to the energy-dependent trans- 
mission coefficient through the barrier, T(E), which it- 
self depends exponentially on the barrier height, <J>, and 
the barrier thickness, t. This can be written as T(E) oc 
exp[-2n($,E)t], with k($,E) = v / 2m e ($ - E)/h, with 
m e being the electron mass. Clearly if the barrier height 
is different for the different spin-species the current po- 
larizations will increase exponentially with the insulating 
layer thickness, leading to full spin polarization for thick 
barrier s 3 ' 12 . Devices constructed with this principles are 
called spin filter tunnel junctions (SFTJs). The quest 
for manufacturing SFTJs then reduces to that of finding 
suitable ferromagnetic semiconductors. 

The europium chalcogcnidcs EuS, EuSc and EuO have 
all a rocksalt structure and they arc all ferromagnetic 
insulators. Among these EuO presents the largest con- 
duction band exchange splitting, ^0.54 eV, below the 
material's Curie temperature of 69 K— ~— . In EuO 
the divalent Eu ions possess a large local moment orig- 
inating from the half- filled 4/ band (//Eu = 7/Ab)- 
An energy gap of 1.1 eV separates the half- filled ma- 
jority Eu-4/ band from the Eu-5g? conduction band 3 -. 
SFTJs based on polycrystalline EuO in the form of a 
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metal/EuO/mctal hctcrojunction have been studied in 
several recent experiments^—. However, the spin trans- 
port properties of crystalline epitaxial EuO junctions 
have not been studied theoretically so far. The pur- 
pose of this paper is to provide such theoretical insight. 
In particular we present ab-initio results for the elec- 
tronic structure and the electron transmission through 
EuO barriers sandwiched between Cu electrodes oriented 
along the [001] direction. 

Our paper is organized as follows. We start our discus- 
sion by presenting the methods used and the structure 
of the device investigated. Then, in the following sec- 
tion we present the electronic structure of EuO and its 
complex band structure along the [001] direction. This 
determines the spin-dependent decay rates and thus the 
spin filter efficiency. In section IV wc discuss the trans- 
mission coefficient at zero bias, while in section V the 
dependence of the spin transport on the EuO thickness 
is analyzed and related to the complex band structure. 
The current-voltage, I-V, characteristics are discussed in 
the following section, where we find that, in contrast to 
Fc/MgO junctions, the polarization of the current is not 
reduced by the applied bias. Finally we conclude. 



II. METHODS 



The EuO and Cu electronic structures are calculated 
by using the ab-initio density functional theory (DFT) 
code SIESTA^. Since the local density approximation 
(LDA) is not sufficient to give an accurate description of 
the EuO density of states (DOS), wc correct for on-site 
Coulomb re puls ion with an LDA+U treatmen t 24 ' 25 . Fol- 
lowing Ref. |26| , the exchange constant J and the on-site 
orbital potential U for the Eu-4/ orbitals are set respec- 
tively to Jf = 0.77 eV and U f = 8.3 eV, while for the 
0-2p orbitals we use J p = 1.2 eV and U p = 4.6 eV. 
Troullier-Martins norm-conversing relativistic pscudopo- 
tcntials are used for Cu, Eu and O. The wave functions 
are expanded over a double £ plus polarization basis set 
(DZP), except for the Eu-4/ states, for which we use 
single-C- An equivalent real space mesh cutoff of 600 Ry 
is used together with an electronic temperature of 69 K. 
We sample the Brillouin zone (BZ) in the plane perpen- 
dicular to the transport direction over a 7 x 7 fc-point 
uniform mesh. 

Spin transport is studied by using the SMEAGOL elec- 
tronic transport code^Zr— , which combines DFT with 
the nonequilibrium Green's functions (NEGF) technique. 
SMEAGOL uses the Hamiltonian matrix provided by 
SIESTA to calculate the charge density so that the same 
pseudopotentials and exchange and correlation function 
can be used for both the electronic structure and the 
transport. The spin current at each bias voltage is calcu- 
lated by energy integration of the spin-dependent trans- 
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where a is the spin index (a =^,^), / is the Fermi 
function, V is the applied bias voltage and e the 
electron charge. The total transmission coefficient is 
obtained by integrating the fc-dependent transmission 
T a (E, V, k x , k y ) over the 2-dimensional BZ perpendicular 
to the transport direction, 



T a (E, V) = YlJ dk x dk y T°{E, V, k x ,k y ) 



(2) 



where SI is the area of the BZ. We denote T(E) = 
T(E, V = 0) as the zero bias transmission coefficient. In 
our calculations wc assume the lattice structure to be pe- 
riodic in the x-y plane, and we keep the transport direc- 
tion along the z-axis. While 7x7 fc-points are enough to 
accurate converge the charge density, the presence of res- 
onances in T' T (E, V, k x , k y ) requires the much larger mesh 
of 100 x 100 to evaluate the transmission coefficient both 
for zero and finite bias. The I-V characteristics are calcu- 
lated non-self-consistcntly by evaluating the transmission 
coefficient over an effective bias-dependent Hamiltonian 
matrix, which in turns is obtained by adding a rigid shift 
to the zero bias Hamiltonian matrix elements of the elec- 
trodes and a linear potential across the insulating barrier. 
This is a good approximation to the self-consistent poten- 
tial drop for tunnel junctions^, which appears essentially 
like that of a standard parallel plate capacitor. 

EuO crystallizes in the rocksalt structure with a lattice 
constant of oeuO = 5.144 A. The primitive face-centered 
cubic (fee) unit cell containing one Eu and one O atom is 
shown in Fig.flja). Since we consider transport along the 
[001] direction, we have to use as basic building block for 
the EuO spacer a cell with size and orientation different 
from that of the primitive one. The smallest possible cell 
thus has tetragonal symmetry, and contains 2 Eu and 2 
O atoms [see Fig. QJb)]. By stacking multiples of such 
tetragonal cells along the z-direction we can construct 
barriers of arbitrary thickness, where one cell contains 2 
EuO monolayers (MLs). 

The model spin filter junction considered here consists 
of n MLs of EuO sandwiched between non-magnetic Cu 
electrodes, n being an integer. As electrodes material we 
consider fee Cu oriented along the [001] direction (lat- 
tice constant 3.61 A), a material choice which has been 
adopted in several experiments^— . When the tetrago- 
nal EuO unit cell is used the lattice constant of Cu can be 
matched with only a slight strain to the EuO lattice, since 
the dimensions along x and y are aEuO / V^2 = 3.63 A. The 
junction setup is illustrated in Fig. [5J In this basic setup 
both the O and Eu atoms are placed over the hollow sites 
of the Cu surface. We also perform calculations for EuO 
shifted in the x-y plane. However the transport and spin 
filter properties arc not sensitively dependent on such a 
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FIG. 1: (Color on line) Lattice structure of EuO constructed 
respectively with (a) the primitive fee unit cell and (b) the 
tetragonal unit cell. The shaded grey atoms indicate the full 
cubic cell. Atom code: large spheres (blue) = Eu, small 
spheres (red)= O. 
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FIG. 2: (Color online) Supercell used for the spin filter tunnel 
junction. This consists of 6 Cu MLs (left lead), 9 Eu MLs 
(scattering region) and 5 Cu MLs (right lead). Atom code: 
large spheres (blue) — Eu, medium spheres (brown) = Cu, 
small spheres (red) = O. 



shift. The equilibrium interface distance, d, between Cu 
and EuO is found to be d = 2.8 A. Note that we do not 
consider possible oxidation at the Cu/EuO interface, or 
the formation of secondary EuO phases. 



III. EuO ELECTRONIC STRUCTURE 

Let us start our discussion by describing the electronic 
structure of EuO. The calculated EuO band structures 
for both the primitive fee and the tetragonal unit cells are 
presented respectively in Figs. [3] and [5] The results agree 
rather well with previously published calculation s 26 ' 30 . 
In order to compare the band structures for two cells 
we use the same /c-space path. Standard high-symmetry 
points are specified for the primitive fee BZ and also for 
the primitive tetragonal BZ (where applicable). We note 
that due to band-folding for the tetragonal cell we find 
twice the number of bands than for the primitive fee 
one. The X point [located at (7r/aE U O,0,0)] and the X' 
point [located at (0,0,7r/aE U o)] are equivalent in the fee 
BZ. However for the tetragonal cell they are independent 
since the X' point is equivalent to T in the tetragonal cell. 
There is a band gap of about 1.0 eV for majority spins 
(Fig. whereas the gap increases to about 3.5 eV for 
the minority ones. For the primitive fee cell (Fig. [3]) we 
find an indirect gap, where the top of the valence band is 
at r and the bottom of the conduction band at X' point. 
This gap become direct when the primitive tetragonal 
cell is used due to band- folding (Fig. [5]). Note that the 
tetragonal cell is the relevant one for the transport so 
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FIG. 3: LDA+U band structure of EuO for the primitive 
fee unit cell [see Fig. [Ha)]: left panel majority spin bands; 
right panel minority spin bands. The band structure shows 
an indirect band gap from f to X of about 1.0 eV. 



that EuO along the [001] direction behaves as a direct 
gap ferromagnetic semiconductor. 

The majority spins the energy spectrum consists of 
three parts: the conduction bands, the top valence bands, 
and the lower valence bands. The top valence bands are 
found in the range between —1.6 eV and —0.5 eV (we set 
the energy zero approximately to the middle of the gap). 
It is clearly seen from Fig. 01 where we show the pro- 
jected DOS (PDOS) on different orbitals, that the main 
character of the top valence band is given by the Eu-4/ 
states. While many of the Eu-4/ bands show very little 
dispersion, indicating strong localization, in EuO the Eu- 
Af/0-2p mixing leads to a significant dispersion around 
the r point. This is an indication of delocalized states. 
We note that at the top of the valence band Eu-4/ and 
0-2p contribute approximately equally to the density of 
states. In addition, the hybridization leads to a large con- 
tribution to the inter-site coupling of kinetic processes 3 -^. 
We therefore expect such delocalized states to contribute 
significantly to the transport through a EuO barrier. Be- 
low the top valence bands, the bottom valence bands are 
observed in the energy range from —4.9 eV to —3.0 eV. 
It is seen from Fig. [3] that the bottom valence bands are 
dominated by 0-2p states. The bandwidth of the 0-2p 
states (1.9 eV) is not much larger than the one of the Eu- 
4/ ones (1.1 eV), which indicates a similar wave- function 
delocalization. The energy gap between the top valence 
and bottom conduction bands is about 1.7 eV. The lower 
conduction bands are dominated by Eu-5g? states. Above 

I. 6 eV we find also contributions from the Eu-6s orbitals. 
We note that the DOS at the lower end of the conduction 
bands is very small, which is due to a high dispersion. 

Next we consider the electronic band structure of the 
minority spins. The empty Eu-4/ states start at around 

II. 1 eV, so that they are not expected to affect the trans- 
port properties. As for the majority states, the conduc- 
tion bands are dominated by Eu-5d and Eu-6s orbitals, 
while the valence bands are dominated by 0-2p. The spin 
filter character of EuO is due to a spin splitting of the 
bands, which leads to a difference in the band gaps for 
the majority and minority spins. The splitting between 
the majority and minority spin conduction bands is cal- 
culated to be about 0.6 eV, in agreement with previous 
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FIG. 4: (Color online) EuO density of states projected over 
the following atomic orbitals: 0-2p, Eu-6s, Eu-5g( and Eu-4/. 
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FIG. 5: LDA+U band structure of EuO for the primitive 
tetragonal unit cell [see Fig. [TJb)], with majority spin bands 
in the left panel and minority spin bands in the right panel. 
Due to the doubled size unit cell, as compared to the primitive 
fee one, a band down-folding makes the T point and the fee 
X' point equivalent. This leads to a direct gap at F. 



theoretical studie d 26 ' 30 " — and with recent experiments 
using 3D angle resolved photoemission spectroscop y 35 ' 36 . 
The Eu-6s states show a similar spin splitting as the Eu- 
5d states, namely 0.4 eV. 

We now analyze the symmetries of the different bands. 
Symmetries are important for the transport through epi- 
taxial tunnel junctions^ 3 -, since, at any given energy, 
only evanescent states in the insulating barrier with 
matching symmetry with propagating Bloch states in the 
electrodes can contribute to the conductance. Since the 
EuO gap in the tetragonal cell is at T, the transport will 
be dominated by states close to k x = k y = 0, denoted as 
Eo,d, for which the barrier height is smallest (the coor- 
dinate system is defined in Fig. [1} . In the energy range 
comprised between —1.8 eV and 1.6 eV around Ep the 
Cu electrodes have only states of Ai symmetry with re- 
spect to the z-axis. For a cubic space group, the Ai 
symmetry transforms as a linear combination of 1 (s- 



orbitals), z (p z orbitals), 2z 



■y (d z 2 orbitals), and 



z(2z 2 - 3x 2 - 3y 2 ) (/ z 3 orbitals). Above 1.6 eV the Cu 
states have A5 symmetry. The A5 symmetry transforms 
as a linear combination of x (p x orbitals), y (p y orbitals), 
xz (d xz orbitals) and yz (d yz orbitals). At energies be- 
low — 1.8 eV we find the Cu-3g? orbitals, so that there 
are also states with different symmetries. For a trans- 
port measurement up to about 2 V we therefore expect 
the states with Ai symmetry to determine the transport 
properties. 

In order to investigate in more details the propagating 
and evanescent states in the EuO barrier we calculate the 
complex band structure (Fig. [5]) at r2D^ ' 13 ' 37 and ana- 
lyze the symmetries of the evanescent states. Real wave 
vectors (Im[/c z ] = 0) represent propagating states and 
complex wave vectors (Im[fc z ] 7^ represent evanescent 
states, since their wave functions decay as exp(— Im[fc z ] -z) 
across the barrier. For the majority spins the top of 
the valence band is three-fold degenerate, with one state 
with Ai symmetry (Eu-4/ z 3 and 0-2p z orbitals) and two 
states with A5 symmetry. Since the state with Ai sym- 
metry has a lower effective mass than those with A5 
symmetry, the corresponding Ai evanescent states have 
smaller Im[fc z ], and therefore a slower decay. The bot- 
tom of the conduction band is given by a state with A2 
symmetry (F,u-5d x 2_ y 2 orbital). At 1.6 eV wc find the 
Eu-4s states, which have Ai symmetry. Therefore Im[/c z ] 
for the Ai band forms a semi-circle between the Eu-4/ z 3 
states at the top of the valence band and the Eu-4s states 
at 1.6 eV. 

Consequently, we expect the tunneling transmission to 
be dominated by the Ai states in this energy range. The 
Eu-5d propagating states with A2 symmetry cannot cou- 
ple to the Ai Cu states, and therefore are not expected to 
contribute significantly to the transmission. For k x and 
k y different from zero the Eu-5g? states are no longer fully 
orthogonal to the Cu-4s states. However, since the gap 
increases with increasing k x and k y , the barrier height for 
such states is larger. Hence their contributions to the to- 
tal transmission are smaller. Overall one expects a rather 
weak coupling to the Cu-4s states and a low transmission 
for states at the bottom of the EuO conduction band. 

In the minority spins bands there are also three de- 
generate states at the top of the valence band, one with 
Ai symmetry (0-2p z orbital) and two with A5 symmetry 
(0-2p x and 0-2p y orbitals). The Ai state again has lower 
effective mass and therefore a smaller decay for evanes- 
cent states. At the bottom of the conduction band the 
symmetries are analogous to those of the majority spins. 
We therefore also expect the minority conductance to be 
dominated by the Ai symmetry. However, since the mi- 
nority band gap is much larger than that of the majority 
spins, Im[fc z ] is also significantly larger for Ai states in 
the gap. As a consequence, we expect to be substan- 
tially larger than T"k This is usually attributed only to 
different energies of the bottom of the conduction band 3 -, 
whereas our complex band structure analysis indicates 
that the decay of the evanescent states has significant 
contributions also from states in the valence band, espe- 
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FIG. 6: (Color online) Complex band structure of EuO in 
its tetragonal cell [see Fig. H^a), where the symmetries are 
indicated for the bands relevant to our transport setup. Red 
lines = minority spins, black lines = majority spins. 



cially for the majority spins. We note that the Eu-5g? A2 
states will contribute to the transport if different elec- 
trodes, possessing A2 symmetry states in the relevant 
energy range, arc used 



IV. SPIN TRANSPORT PROPERTIES OF THE 
EuO JUNCTION AT ZERO BIAS 

In this section wc analyze the zero bias transport prop- 
erties of Cu/EuO/Cu junctions by taking as an example a 
stack containing 9 MLs of EuO (t = 20.58 A). The energy 
level alignment between the metal and the insulator is an 
important factor determining the transport properties of 
the junction. To a first approximation, one can estimate 
the alignment by comparing the workf unctions, W, of the 
two components. The workfunctions of EuO and Cu are 
calculated by using the Hartree electrostatic potential, 
Vh, as the reference potential^. We define Wc u as the 
difference between Ep and the vacuum potential Vv acuum 
of a Cu slab, while WeuO is given by the difference be- 
tween Vvacuum and the energy of the valence band top 
of a EuO slab, Vvb- We calculate Wc u = 3.9 eV and 
WeuO = 1-8 eV, both values in good agreement with the 
experimental values, respectively ranging between 4.5 eV 
to 5.1 eV for Cu, and being 1.7 eV for Eu02£ Since 
W-EuO < Wcu, electron transfer from EuO to Cu can be 
expected at the interface, leading to the pinning of the 
valence band top of EuO to the Cu Ep. This is indeed 
the case for the equilibrium distance (d = 2.8 A), where 
the Eu-4/ states are located just below Ep [see Fig. [TJc) 
and (d)]. 

In experiments investigating bias-dependent transport 
the actual junctions usually contain polycrystalline EuO, 
so that the interface between Cu and EuO is not well 
defined 22 . In our calculations, in contrast, we assume 
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FIG. 7: (Color online) PDOS for Cu, O and Eu for a junction 
comprising 9 MLs EuO. The PDOS is calculated at the center 
of the junction (a) and the interface (b) for d = 2.4 A, and 
at the center of the junction (c) and the interface (d) for 
d — 2.8 A. Here d is the distance between the Cu and the 
EuO planes at the EuO/Cu interface. 



a perfect epitaxial interface and do not explicitly con- 
sider the formation of defects, oxidation of Cu, or the 
possible formation of EU2O3 at the interfac o 17 ' 21 . The 
effect of such modifications of the interface on the elec- 
tronic structure are manyfold, the most important being 
that they usually lead to a different charge transfer and 
therefore to a different band alignment between Cu and 
EuO. In practice this means that (depending on the de- 
tailed structure of the interface, which is determined by 
the experimental conditions) Ep can be placed at differ- 
ent positions across the EuO gap. In our calculations we 
can tune the charge transfer and consequently the posi- 
tion of Ep in the EuO gap by modifying the distance, 
d, between Cu and EuO at the interface. We find that 
decreasing the distance from the equilibrium one, i.e., 
increasing the coupling between Cu and EuO, the EuO 
states shift to lower energies with respect to Ep. For 
d = 2.4 A Ep is located approximately in the middle 
of the gap [see Fig. [TJa) and (b)], while for d = 2.2 A 
it is pinned at the EuO conduction band minimum (not 
shown). Therefore, d = 2.2 A can be used to simulate 
the transport properties of n-type EuO, as obtained for 
0-dcficicnt barriers^—. In the remaining part of this 
paper we will present the transport properties for both 
d = 2.4 A and d = 2.8 A, in order to illustrate the ef- 
fect of a shift of Ep induced by interface modifications. 
Importantly, we will demonstrate that, for any position 
of Ep, EuO always shows excellent spin filter character- 
istics up to high bias voltages. We note that if experi- 
mentally a perfect epitaxial junction can be realized, the 
measurements should correspond to our results for the 
equilibrium distance d = 2.8 A. 

In Fig. [7] we compare the PDOS of atoms at the 
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FIG. 8: Transmission coefficient of a Cu/EuO/Cu junction 
formed by 9 EuO MLs for d = 2.4 A and d = 2.8 A, plotted 
on a linear (a) and on a logarithmic scale (b). Note that for 
d = 2.8 A, T(E) shifts to higher energy when compared to 
the case of d = 2.4 A. 



Cu/EuO interface with the PDOS of atoms at a maximal 
distance from the interface (to be considered bulk- like). 
While the two PDOSs arc similar for most atomic or- 
bitals, there is a significant difference for the 0-2p states. 
At the interface the 0-2p states extend over a larger en- 
ergy range as compared to bulk-like atoms, indicating a 
coupling to the Cu substrate. In fact, for d = 2.4 A the 
broadening is more pronounced than for d = 2.8 A, which 
is due to the larger coupling. The PDOS in the middle 
of the EuO layer is very similar to that of bulk EuO, in- 
dicating that in the middle of the junction one recovers 
the bulk electronic structure of EuO. Although in Fig. [7J 
only the fully occupied Cu-3d states are shown, we note 
that the Cu-4s states determine the transport properties 
of the electrodes, since they have an approximately con- 
stant PDOS in a large energy range around Ep. In Fig. 
[7] the Cu-4s PDOS is not shown, since it is not visible on 
the chosen scale. 

The zero bias transmission coefficient T(E) for differ- 
ent d is shown in Fig. [5] [on a linear scale in panel (a) 
and on a logarithmic scale in panel (b)]. Due to the shift 
of the electronic states for different d (see Fig. UJ), T(E) 
is shifted towards higher energies for d = 2.8 A as com- 
pared to d = 2.4 A, while there is only a minor change in 
the height of the different transmission peaks. It can be 
seen that T(E) has two small gaps for the majority spins 
and one large gap for the the minority ones. For the 
minority spins the conduction occurs only through the 
0-2p states for energies below —2.5 eV and through the 
Eu-5d and Eu-6s states for energies above about 2 eV. In 
contrast, for the majority spins the Eu-4/ states just be- 
low Ep also contribute significantly to the conductance. 
This is consistent with the band structure, which shows 
that the Eu-4/ states hybridize with the 0-2p and are 



rather delocalizcd. The result of such hybridization is a 
large majority transmission in the energy range between 
— 1.5 eV and Ep. Therefore, for all energies below Ep 
(down to about —1.5 eV) the transmission of the major- 
ity spins is much larger than that of the minority ones. 
For energies above Ep (up to about 2 eV) it is also sig- 
nificantly larger. This is due to fact that the Eu-5d con- 
duction band minimum is located about 0.6 eV lower in 
energy for the majority spins than for the minority ones 
(Fig.®. 

The overall result hence is that for any position of Ep 
in the EuO gap we expect a very high spin-polarization of 
the current. This should persist up to high bias voltages 
(of the order of the energy for which the transmission of 
the majority spins is much larger than that of the minor- 
ity ones). The spin- filtering efficiency at small applied 
bias is defined as [T*(E F ) - T^(E F )}/T{E F )&. Since 
T^{E F ) is about two to three orders of magnitude larger 
than T^(E F ), we have a spin filter efficiency close to 
100%, indicating that the EuO barrier is an almost per- 
fect spin filter. We note that the extremely high efficiency 
is obtained for defect-free epitaxial junctions. Such value 
might be reduced by defects in the EuO barrier as well as 
for polycrystalline EuO. Nevertheless one can expect that 
also in these cases the spin-filter efficiency may remain 
high. 

In order to investigate the dependence of the spin 
transport properties on the geometry of the Cu/EuO in- 
terface we calculate T(E) for different positions of EuO 
with respect to Cu. So far we have considered a setup 
in which O is placed on top of the hollow site of the Cu 
surface. We compare the results with the following three 
geometries: 1) O placed on top of Cu, 2) O placed on the 
bridge site between two Cu atoms, and 3) O placed at an 
arbitrary low symmetry site. The equilibrium distance 
is 2.8 A for all cases and the total energies are similar, 
with the bridge site having the lowest energy. In Fig. O 
we show T(E) calculated for the different sites, for both 
d = 2.4 A and d = 2.8 A. Overall the changes in transmis- 
sion corresponding to the different sites are rather small, 
a fact, which indicates that the geometry of the Cu/EuO 
surface does not significantly affect the band alignment 
and transport properties of the junction. 

At the bottom of the conduction band there are some 
quantitative differences in transmission between the dif- 
ferent sites. This is due to the fact that shifting the 
O atom with respect to Cu alters the coupling between 
the Cu-4s and Eu-5d states. As discussed in Section Hill 
while the states on the top of the valence band couple 
well to the Ai Cu band, the states at the bottom of 
the conduction band couple only slightly to these states 
due to a symmetry mismatch. This is also the reason 
for the fact that the transmission gap is somewhat larger 
than the EuO band gap, and for the slow increase of 
the transmission at the bottom of the conduction band 
with energy. For the bridge and low symmetry sites the 
symmetry mismatch is slightly reduced, resulting in a 
somewhat smaller transmission gap. Since the results for 
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FIG. 9: (Color online) Transmission coefficient for a 
Cu/EuO/Cu junction with a 9 MLs EuO barrier and with 
different interface geometries: (a) d = 2.4 A and (b) d = 2.8 
A. The solid lines represent the majority spin transmission, 
whereas the dashed ones are for minority spins. Note that 
overall the transmission depends little on the lateral position 
of the O atoms with respect to the Cu surface. 



the different sites are very similar, we will only consider 
the hollow site in the remaining of the paper. 



V. THICKNESS DEPENDENCE OF THE 
CONDUCTANCE 

In this section we analyze the effect of the EuO 
thickness, t, on the spin transport properties of the 
Cu/EuO/Cu junction. The transmission spectra for dif- 
ferent EuO thicknesses, and for both d — 2.4 A and 
d = 2.8 A, are shown in Fig. [TO] From the figure it can 
be seen that the band alignment is not affected by the 
thickness, while the transmission in the gap decreases ex- 
ponentially with it. Since for energies in the gap we have 
T a (E) ex exp[— 2K a (E)t], we can estimate k ct from the 
change of T a (E) with increasing t. Indeed kF can be cal- 
culated from T a evaluated at two different thicknesses t\ 
and t2 as 
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2{t 2 -t x ) 



In 



THE) 
T%{E) 



(3) 



where T" and T-f are the transmissions for t\ and t 2l 
respectively 

In Fig. [IT] we show the calculated k(E) (solid lines), 
evaluated for t x = 20.58 A (9 EuO MLs) and t 2 = 25.72 A 
(11 EuO MLs). For energies within ± 1.5 eV around Ef, 
k is significantly smaller than ktK Therefore, increasing 
the EuO layer thickness leads to an enhancement of the 
ratio between the transmissions of the majority and mi- 
nority spins, i.e., to an increase of the spin filtering effi- 
ciency For d ~ 2.4 A, is larger than zero in the range 
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FIG. 10: (Color online) Thickness dependence of T(E) at zero 
bias for (a) d = 2.4 A and (b) d = 2.8 A. The solid and dashed 
lines refer to the majority and minority spins, respectively. At 
d = 2.4 A the tunneling gap of the majority spins is found 
from —0.45 eV to 1.1 eV, while for the minority ones it extends 
from —2.5 eV to 1.6 eV. At d = 2.8 A the tunneling gaps of the 
majority and minority spins for all thicknesses range between 
-0.1 eV to 1.75 eV and -2.3 eV to 2.5 eV, respectively. 
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FIG. 11: Damping coefficient n calculated from T(E) for the 
majority [(a) and (b)] and minority spins [(c) and (d)] at 
d = 2.4 A and d = 2.8 A. k is calculated from Eq. [3] by 
using the transmission coefficient calculated for 9 and 11 EuO 
MLs. The solid lines refer to the total transmissions and the 
dashed lines to the transmissions at k x = k y — only. The 
green lines show Im[fc z ] from Fig. [6] which is also calculated 
for k x — k y = 0. 



from —0.5 eV to 1.0 eV (resulting in a transmission gap 
of 1.5 eV), whereas is larger than zero in the range 
from —2.6 eV to 1.6 eV (resulting in a transmission gap 
of 4.2 eV). For d= 2.8 A, k is shifted to higher energies, 
and the transmission gap amounts to about 1.7 eV for the 
majority spins and 4.7 eV for the minority ones. While 
k shows a parabolic behavior for energies close to the 
valence band top, for the conduction band minimum the 
behavior is less well defined. This is due to the symmetry 
mismatch between the Cu Ai states and the conduction 
band Eu-5d states, see Sec. IIVI Nevertheless the barrier 
height is about 0.6 eV larger for the minority spins as 
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compared to the majority ones for transmission through 
the conduction band. For the transmission through the 
valence band the difference is even larger, since there are 
no filled minority spin Eu-4/ states. 

In Fig. [11] we also present the k(E) obtained from the 
transmission only at the point (dashed curve). Close 
to the valence band the wave-function decay is similar to 
the previous case, which shows that T(E) is dominated 
by contributions around T2D- However, for energies close 
to the conduction band, k(E) at is much larger than 
the total k(E), which shows that here the transport oc- 
curs mainly at fc-points away from the BZ center. A 
comparison with the EuO complex bands (green curves; 
see Fig.[B]for a description of the symmetries) shows that, 
at T2D, K (E) follows approximately Im[fc z ] for Ai states. 
It is clear that no transmission occurs through the con- 
duction band A2 states. For fc-points away from T2D the 
EuO A2 states can couple to the Cu states, which leads 
to a decrease of the total k(E) above the conduction band 
minimum. For electrode materials with A2 states above 
Ep we expect the total k(E) to follow the A2 complex 
bands at I^d for energies below the conduction band 
minimum. 

The barrier heights are calculated from the damping 
coefficient shown in Fig. [TT] The values of and 
are determined as the lowest energy at which k(E) crosses 
zero in the conduction band for the majority and mi- 
nority spins, respectively. Similarly, $y B and <I>y B are 
determined as the highest energy at which n(E) crosses 
zero in valence band for the majority and minority spins, 
respectively. By this method we obtain for d = 2.4 A 
the values = 1.0 eV, <f>^ B = 1.6 cV, $ VB = 0.5 cV, 
and <f>^ B = 2.6 eV. For d = 2.8 A we obtain $£ B = 1.6 



cV, <f>£ B = 2.2 eV, ${t3 = 0- 1 cV > and $ vb = 2 -4 eV. 
The exchange splitting for the conduction band of 0.6 eV 
matches the value for bulk EuO. 

For the low bias conductance the thickness dependence 
of the transmission coefficient at Ep is evaluated and 
from it the decay coefficient. In Fig. [12] we show T(Ep) 
at d = 2.4 A and d = 2.8 A as functions of the EuO 
thickness. The decay rates correspond to those obtained 
in Fig. [TT] at Ep. In Fig. [T2Tc) the effective spin filter- 
ing at Ep is addressed. It can be seen that it increases 
towards 100% as t increases. 



VI. THE I-V CURVE OF THE EuO SPIN 
FILTER 
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FIG. 12: Dependence of T(E F ) on t for (a) d = 2.4 A and (b) 
d = 2.8 A. The spin filtering efficiency at Ef is shown in (c) 
for both values of d. 
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The spin-resolved current- voltage, I-V , characteristic 
and the spin polarization as functions of the bias voltage 
are shown in Fig. Q2] for 7, 9 and 11 MLs of EuO, and 
for both d = 2.4 and 2.8 A. The majority spin current is 
higher than the minority spin current for all bias voltages 
and for all thicknesses, leading to large polarization in all 
cases. At low bias there is a rapid increase of the tunnel- 
ing current with the bias. This is typical for such systems 



FIG. 13: Spin-polarized current and spin polarization as func- 
tions of the bias voltage for different EuO thicknesses, (al) 
Majority spin current at d = 2.4 A, (a2) majority spin current 
at d — 2.8 A, (bl) minority spin current at d = 2.4 A, (b2) 
minority spin current at d = 2.8 A, (cl) spin polarization at 
d — 2.4 A, and (c2) spin polarization at d = 2.8 A. 
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FIG. 14: Total current and conductance for a Cu/EuO/Cu 
junction with a 9 EuO MLs barrier at d = 2.4 and 2.8 A with 
different interface geometries, (al) Current at d = 2.4 A, 
(a2) current at d = 2.8 A, (bl) conductance at 2.4 A, and 
(b2) conductance at 2.8 A. 

due to the fact that the current is not only determined 
by electrons at Ep but also by those in the energy range 
Ep ± eV/2 (bias window), for which the barrier height 
can be reduced up to about eV/2. At a bias of about 
0.7 V for d = 2.4 A the electrons start to flow through 
the valence band top, so that there is a sharp increase 
in current. Since for d = 2.8 A the valence band is very 
close to Ep, the onset is found at very low bias. The 
current increases up to about 1.3 V, above which a de- 
crease with increasing bias is found. This is caused by 
a reduction of the transmission over the entire bias win- 
dow once the bias is very large. Such a current reduction 
results from the enhanced scattering as the potential is 
tilted inside the barrier due to the applied bias voltage. 
In these large scattering conditions, incoherent contribu- 
tions to the current (not included here) are expected to 
play an important role. 

For the minority spins, in contrast, the current remains 
in the tunneling regime for all bias voltages due to the 
large barrier height. The large difference between major- 
ity and minority spin currents is reflected by a high spin 
polarization, P, defined as (I* — 7^) / (U + 1^) [see Fig. 
[T3T cl)]. P remains large for all the bias voltages consid- 
ered due to the small 1^. We note that the spin splitting 
of the conduction band plays only a secondary role in our 
results, since the majority current is mainly determined 



by the valence band contribution (which is not sensitive 
to the exact position of Ep inside the EuO energy gap). 
For n-typc EuO one expects a pinning of the conduction 
band to Ep, leading to an almost metallic-like character 
dominated by the EuO conduction electrons. 

The dependence of the I-V characteristics and of the 
conductance trace G(V) = dl/dV on the Cu/EuO inter- 
face geometry is addressed in Fig. [TU Overall the I-V 
characteristics are similar for all the interfaces, which 
reflects the fact that the zero bias transmission is also 
rather similar (see Fig. [9]). For d = 2.4 A there is a sharp 
increase in G at the current onset, whereas it then stays 
approximately constant. For d = 2.8 A, in contrast, G is 
large at low bias and changes sign at 1.3 V. 

VII. CONCLUSIONS 

We have studied the electronic structure, the complex 
band structure and the spin transport properties of epi- 
taxial Cu/EuO/Cu tunnel junctions. The spin transport 
properties of EuO epitaxially grown on Cu are dominated 
by the Eu-4/ valence states and by the Eu-5d conduc- 
tion states. We show that EuO acts as an almost perfect 
spin filter, where close to 100% spin polarization can be 
achieved. The polarization increases with increasing EuO 
thickness as expected from the complex band structure, 
where the decay of the wave function into EuO is pre- 
dicted to be much smaller for the majority spins than 
for the minority ones. Since the conduction states of 
EuO have no Ai symmetry to match the states in the 
Cu electrodes, the gap in the transmission coefficient is 
significantly larger than the band gap of EuO. Under a 
bias voltage the spin polarization of the current does not 
decrease, as in usual tunnel junctions, but remains ap- 
proximately constant up to all considered bias voltages. 
This opens promising perspectives of using EuO for de- 
vice applications. 
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